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ABSTRACT 

Using star counts in both color and magnitude from the Two Micron All Sky Survey (2MASS) Second 
Incremental Release Point Source Catalog we search for evidence of non-uniform extinction and stellar 
population density changes in the Galactic Plane. Extinction causes the entire main sequence to shift 
toward redder colors on a color magnitude diagram. A local increase in the stellar density causes an 
increase in the star counts along a line parallel to the main sequence. We find streaks in star count 
color magnitude contour plots along the angle of the main sequence which are likely to be caused by 
distant gas clouds and stellar density variations. The distance of a gas cloud or stellar density change 
can be estimated from the location of the shift in the star count contours. We identify features in these 
diagrams which are coherent across at least 10 degrees in Galactic longitude. A series of features is seen 
at the plausible distance of the expected Perseus spiral arm at a distance of 2 to 4 kpc from the sun. 
However other features as prominent are found at both at larger and smaller distances. These structures 
are over 300 pc in size and so likely to be associated with large scale coherent structures in the gas 
distribution such as weak spiral arms. The presence of multiple and weak spiral arms, and lack of strong 
ones suggests that the outer Milky Way is flocculent in its morphology. 

Subject headings: starcounts 2MASS spiral structure Milky Way 



1. INTRODUCTION 

It has long been recognized that star counts or the num- 
ber of stars per magnitude bin per unit area on the sky are 
strongly affected by extinction, or the absorption of light 
by dust (e.g., Seeliger 1900; McCuskey 1965). Because ex- 
tinction is particularly strong in the plane of the Milky 
Way, it has not been possible to study the spiral struc- 
ture of our Galaxy deep within the Galactic Plane using 
star counts. However, above the scale height of the gas 
and dust, at Galactic latitudes outside the "Zone of avoid- 
ance" or greater than 10° from the Galactic Plane, it is 
possible to measure the distribution of stars in the galaxy. 
For example star counts from the Sloan Digital Sky Sur- 
vey well above the Galactic Plane have recently been used 
to measure the scale heights of the thin and thick stellar 
disks (Chen et al. 2001) and uncover structure in the halo 
(Yanny et al. 2000). 

Because of the sensitivity to extinction, rather than 
probe for galactic structure, star counts are often used 
to measure extinction, as first illustrated with the Wolf 
diagram. Wolf 1923; Trumpler & Weaver 1953, or make 
extinction maps in the visible bands (e.g., Dickman 1978) 
and at very high levels of extinction in molecular clouds in 
the near-infrared (e.g., Lada et al. 1994; Cambresy 1999; 
Lombardi & Alves 2001). 

The 2MASS team has recently released a significant 
fraction of the sky at a depth of about 15 mag in Kg band 
(the Second Incremental Release). For information on the 
2MASS Second Incremental Release Point Source Catalog 
see http://irsa.ipac.caltech.edu. Because 2MASS is car- 
ried out between 1 and 2 microns, the effect of extinction 
is greatly reduced compared to that at visible wavelengths. 
2MASS therefore provides us with a unique opportunity 
to probe for structure within the plane of our Galaxy. 



2. NUMBER COUNTS IN A COLOR-MAGNITUDE DIAGRAM 

To detect structure in the halo Yanny et al. (2000) se- 
lected stars in a narrow range of color. On the main se- 
quence color effectively determines the luminosity of the 
star and so its distance can be determined from the ob- 
served magnitude. While this technique worked well to un- 
cover structure in the halo, we now explain why it presents 
difficulties when used to search for structure in the Galac- 
tic Plane using 2MASS data. 

A difference in J — of 0.1 results in a change in abso- 
lute magnitude on the main sequence of about a magnitude 
(see Fig. 1) and this would correspond to a factor in dis- 
tance of 1.6. The nearest spiral arms are about 2 kpc away 
from the sun (e.g., Vallee 1995) and may only have thick- 
nesses of a few hundred pes, so to see an enhancement in 
the star counts we can only tolerate distance uncertainties 
that are less than factors of 1.1. This corresponds to a 
magnitude uncertainty of less than 0.2 which is difficult to 
achieve using an infrared color section for stars at the dis- 
tance where we expect spiral arms. We expect that mag- 
nitude uncertainties resulting from the strong dependence 
of absolute magnitude on color of the main sequence will 
wipe out any evidence for spiral structure in a given color 
selected sample. Furthermore because of reddening asso- 
ciated with extinction, a color selection will not restrict 
the stars to a give spectral type or absolute magnitude. 
So despite the fact that this approach was successful for 
detecting structure in the halo, if we only search in one 
color bin it is unlikely to work in the Milky Way disk us- 
ing 2MASS data. 

So how do we get around this problem? We must use 
information from stars which cover a range of observed col- 
ors. The main sequence provides us with a useful relation 
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between color and absolute magnitude. A change in stellar 
population caused by recent star formation causes a local 
increase in the number of blue stars, whereas extinction 
will cause the entire main sequence to shift toward redder 
colors on a color magnitude diagram. If wc sec features on 
a color magnitude number count diagram (number counts 
per unit color and magnitude bin) that extend along the 
direction of the main sequence, then we expect that those 
features are caused by local changes in either extinction or 
the local stellar population density. 

In Figure 1 we show star counts from a 2.5 square degree 
field taken from the 2MASS Second Incremental Release 
Point Source Catalog. Stars were counted in an array of 
bins defined by both color and magnitude. Magnitude bins 
wcirc! 0.1 mag wide in either J or Ks bands and J — H color 
bins were 0.025 mag wide. To aid the contouring algorithm 
we slightly smoothed the array with a smoothing function 
3 pixels high and wide that is shaped like a pyramid with 
half the flux in the central pixel. The minimum contour 
shown contains 5 stars per bin and wc show a total of 12 
contours spaced logarithmically with the maximum being 
the maximum value of counts in the array. For the 2MASS 
Point Source Catalog the Level 1 requirements were to 
achieve a signal to noise of S/N= 10 at 15.8, 15.1, and 
14.3 mag in J, H and Kg bands respectively. The actual 
achieved limiting magnitudes at this signal to noise are 
about 0.5 magnitude better than expected. Typical pho- 
tometric errors are about 0.03 mag for 13,12.5 and 12.5 
magnitudes in J, H and Kg band respectively and about 
0.1 mag for 16.3,15.5 and 15.0 mag respectively. 

As a stellar population becomes more distant, the lu- 
minosity function is pushed to higher magnitudes on this 
figiir{\ However extinction both reddens and decreases the 
luminosity and so affects the entire main sequence along 
the same vector in the diagram. A region containing a 
young population increases the number of blue stars at a 
particular magnitude range or distance. An increase in 
the stellar density at a particular ciistancic will cause an 
increase in counts again along an angle parallel to that of 
the main sequence. This will cause a shift in the contours 
toward the left. In Figure 1 we have over- plotted the main 
sequence at a solar metallicity and distance of Ikpc for 3 
different aged populations. Diagonal streaks in the con- 
tour levels are observed which are along the same angle as 
the main sequence but offset in magnitude. Because we 
sec bends in the contours that extend over a magnitude 
in J and K and 0.2 mag in color, they are likely to be 
real and not to due Poisson noise. These diagonal streaks 
are probably due to extinction at a distance exceeding the 
main sequence model points. The distance to the gas cloud 
responsible for the extinction can be estimated from the 
magnitude shift of the main sequence. A streak in the 
contours is seen in both J and Kg vs J — H band color 
magnitude plots corresponding to an extended enhanced 
region of extinction or gas cloud at a distance of about 2.5 
kpc from the sun. 

3. PROSPECTING FOR SPIRAL STRUCTURE 

We now look for coherent structures seen in the color 
magnitude diagrams in different locations in the Galactic 
Plane. There are three regions in the Second Incremen- 
tal release that contain the Galactic Plane and not the 



Galactic Center. These regions have Galactic longitudes 
which range glon = 220 250°, 50 80° and 170-190°. All 
three regions should contain the the Perseus Spiral arm 
(see Fig. 2), however in the 50-80° region the arm is dis- 
tant at 5-7 kpc and so it should be impossible to detect it 
using main sequence stars. Furthermore because the spiral 
arms increase with distance with increasing Galactic lon- 
gitude, even weak spiral arms in the region between the 
Carina and Perseus arms at glon ~ 80° should not extend 
over a large range in galactic longitude. In the Galactic 
anti-center region between glon = 170 — 190° the Perseus 
spiral arm should be roughly equidistant from the sun at 
a distance of 2 kpc but little extinction is evident in this 
direction. So we begin by concentrating on the region be- 
tween 220 and 250°. 

In Figure 3 we show color magnitude number count di- 
agrams in regions ranging from 220-250° with Galactic 
latitudes glat = to —2 or —3 to —1°. In a few regions 
pieces of the sky were not covered by the data release so 
the counts were corrected by dividing by the percent of 
sky area covered. The Kg band plot corresponding to the 
J-band plot in Fig. 3a is not shown since it is quite similar 
to Fig. 3a. In this region between we see streaks (de- 
noted by a dotted line in Fig. 2) that begin at J ^ 14.5 at 
glon = 220 and end at J ~ 15.5 at glon ~ 235. This fea- 
ture is also seen in the Kg color mag contour plot (Fig. 3c). 
The location of the streaks is at somewhat larger magni- 
tudes for glat = —3 to —1 than for glat = 2° suggesting 
that there is a larger smooth extinction gradient with ra- 
dius at this Galactic latitude than at glat = 1°. From the 
location of the feature in the ii's-band plots and how it is 
offset from the main sequence at 1 kpc (see Figure 1) we 
can estimate its distance to be ~ 2.5 kpc at glon = 220 
extending to ~ 4 kpc at glon = 235. This corresponds to 
the expected location of the Perseus spiral arm (see Fig. 2) 
as extended from our knowledge in existing spiral tracers 
(Vallee 1995). 

We also see evidence for additional structures which 
have not been identified by Drimmel & Spergel (2000) 
and HI surveys. See for example the indent in contours 
at glon ~ 225 to 235 at Kg ^ 13.8 to 15 which would cor- 
respond to extinction at a distance of 2-4 kpc but inside 
the Perseus arm which is expected to be at a distance of 
5 kpc at glon = 235°. 

Structure seen in these star count color magnitude con- 
tour plots appears to be is correlated and extended. A 
streak in the contours at one galactic longitude is likely to 
be present in the diagram at lower and higher longitude, 
though systematically shifted toward more higher magni- 
tudes at higher galactic longitudes. This makes sense if 
the density perturbations both in the gas and stars are 
spiral in shaped, or increasing with distance with increas- 
ing Galactic longitude (see Fig. 2,4). 

The structure see in these star count diagrams cannot 
be solely due to extinction because local increases in the 
star counts at higher magnitudes is also seen. This im- 
plies that there are also density variations as a function of 
distance from the sun. Again this is not unexpected if we 
remember what other galaxies look like (e.g., see Fig. 4 
which displays a B-band image of MIOO.) 

Though see strong evidence for local density and extinc- 
tion variations in the region around glon ~ 230, we do not 
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features associated with extremely strong spiral structure, 
neither evidence for large variations in the stellar density 
or local extinction. The most likely explanation for the 
structure that we do see is that the Milky Way is floccu- 
lent in its morphology outside the solar circle. 

For a spiral arm outside the solar circle we expect to en- 
counter first an increase in extinction and then an increase 
in the stellar density as a function of distance from the sun 
as a spiral arm is crossed (see Fig. 4). The affect of the 
extinction is therefore decreased by the increased density 
following the spiral arm. Stellar density and extinction 
variations associated with spiral structure should be eas- 
ier to uncover using star counts in the Carina arm within 
the solar circle because we expect that blue stars should be 
nearer than the belt of extinc;tion associated with the spi- 
ral arm. This should emphasize features seen in the color 
magnitude star count plots, rather than decrease them as 
is true when we look toward larger Galactocentric radii. 
Features associated with these spiral arms may also be 
easier to detect because these spiral features should have 
higher gas densities, hence heavier extinction, and stronger 
stellar density contrasts. 

In Figure 5 we show color magnitude star count plots 
taken from a strip along glat — 2-3° and between glon = 
65—85°. There is a lot of extinction at glon ^ 8° which 
corresponds to the Dark Nebula LDN-906 (Lynds 1962) 
which has an area of 15° and must be quite nearby. A 
number of features are seen in these plots which are prob- 
ably real because they extend over a fair range on an in- 
dividual plot. However we do not find a high level of co- 
herence between the plots. This is not surprising since at 
this Galactic longitude we are looking between the Carina 
and Perseus arms and the opening angle of the arms is 
such that elongated gas clouds will not extend over a large 
range in Galactic longitude. The large number of features 
seen in these plots is consistent with a patchy ISM and 
the possibility that the outer Milky Way is flocculent in 
its morphology. 

In Figure 6 we show color magnitude star count plots 
taken from strips along glat = —1 to 1°, —3 to —1° and 
between glon = 170-190°. At glon - 170 to 180° we 
see evidence for extinction at a distance of about 2 kpc 
(streaks beginning at Ks ~ 13.8). This does corresponds 
to the expected location of the Perseus spiral arm, however 
the lack of any features at glon ~ 190 suggests that it is 
not continuous. There are also extended features further 
out at Ks ^ 14.5 from glon ~ 180 to 190°. We see no dom- 
inated feature at expected location of the Perseus arm but 
instead a variety of shorter features at both shorter and 
longer distances. 

4. SUMMARY AND DISCUSSION 

In this paper we have presented color magnitude number 

count contour plots constructed from the 2MASS Point 
Source Catalog in the Galactic Plane. We demonstrate 
how streaks observed in these diagrams c;an be used to es- 
timate the distance of local regions of high extinction. A 
coherent structure is seen between Galactic longitude 220 



and 237° which corresponds to the probably location of 
the extension of the Perseus spiral arm. Another coher- 
ent structure is seen with a steeper pitch angle between 
Galactic longitude 230-240° at a distance of 2-4.0 kpc and 
may be a weaker spiral feature. Neither feature is particu- 
larly prominent. In the region glon ~ 170 to 190° we see a 
feature at the probably location of the Perseus spiral arm 
but also also other coherent features at both larger and 
smaller distances. Features seen in the star count contour 
plots are likely to be real because they cover a significant 
range in color and magnitude in each color magnitude star 
count plot and because nearby lines of sight in the galaxy 
exhibit similar features. At a distance of 2kpc a structure 
which extends over 10 degrees has a length of 350 pc so 
the; coherent structures picked out from the star counts are 
likely to associated with spiral arms. The lack of strong 
spiral features and the evidence for additional weak fea- 
tures suggests that the outer part of Galaxy is flocculent 
in its morphology. 

We suggest that color magnitude number counts could 
be used to estimate the distances to large scale gas clouds 
at distances up to a few kpc using blue main sequence 
stars. It is possible that a similar technique could also 
reveal structure out to much larger distances using gi- 
ant stars. Giant stars have the disadvantage that it is 
not straightforward to estimate the luminosity of the star 
given its color but the advantage that they can be seen 
much further away than blue main sequence stars. 

Using star counts to measure distances to atomic and 
molecular clouds is potentially a very powerful technique. 
However to do this reliably we would need to integrate the 
stellar luminosity function as a function of distance from 
us and devise techniques to differentiate between patchy 
extinction, stellar density variations, mean extinction vari- 
ations and variations in the mean age of the stellar pop- 
ulation. All of these affects should affect the integrated 
number counts. By combining information available in 
CO, HI and far-IR images of the Milky Way with 2MASS 
it may be possible to make a map of the Milky Way with 
the detail of that shown in Figure 4 of Ml 00 rather than 
that currently known which is represented in Figure 2. 

This work could not have been carried out without 
helpful discussions with Larry Heifer, Rob Selkowitz, Joel 

Green, Dan Watson, and Eric Blackman. 

This research has made use of the NASA/IPAC In- 
frared Science Archive, and the NASA/IPAC Extragalac- 
tic Database (NED) which are operated by the Jet Propul- 
sion Laboratory, California Institute of Technology, under 
contract with the National Aeronautics and Space Ad- 
ministration. This publication makes use of data prod- 
ucts from the Two Micron All Sky Survey, which is a 
joint project of the University of Massachusetts and the 
Infrared Processing and Analysis Center/California Insti- 
tute of Technology, fmided by the National Aeronautics 
and Space Administration and the National Science Foun- 
dation. 
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Fig. 1. — a) Color magnitude number counts shown for point sources from the 2MASS point source 
catalog with 222.5 < glon < 225.0 and —3 < glat < — 1° for the color J — H (x-axis) versus J(mag) 
(y-axis). Contours are shown at 12 logarithmic levels with the lowest level set at 10 stars per bin 
and the highest at the maximum of the counts. As data points we show the main sequence at 3 
different ages, 10^, 10^ and 10^'^ years and with solar metallicity assuming a distance of 1 kpc. To 
plot these points we have used the stellar evolutionary tracks of Girardi et al.(2000). An extinction 
arrow is drawn for = 0.5 assuming an extinction law from Mat}iis(1990). As a stellar population 
becomes more distant, the luminosity function is pushed to higher magnitudes. However extinction 
both reddens and decreases the luminosity. A region containing a young population will increase 
the number of blue stars at a particular magnitude range. Diagonal streaks in the contours can 
be caused by extinction or an increase in the stellar number density population at a particular 
distance. Note for example the streak which begins at J « 14.5 and J — H 0.20 and extends 
to J « 15.5 and J — H ^ 0.4. This streak is also seen in b) beginning at K « 14.2 and ending 
Sit K 15. This streak would correspond to enhanced extinction at distance of about 2.5 kpc. b) 
Same as a) except the y-axis is Kg (mag). 



Fig. 2. — Model for the Milky Way spiral arms constructed from dust emission seen in 
COBE/DIRBE data by Drimmel & Spergel(2001). We have overlayed Galactic longitudes at the 
position of the sun so that it is clear what spiral arms we expect to see at a given position on the 
sky and distance from the sun. Shown as a faint gray line passing through the sun is the Orion 
spur, a local ridge of gas and star formation. 




(5duj) p (&OLU) ^yi (boiu) p 



Fig. 3. — a) Color magnitude number count diagrams are shown for J(mag) versus J — H color for 
regions 2.5° wide in galactic longitude and for latitude —3° < glat < — 1°. The leftmost panel has 
220° < glon < 222.5° and the rightmost panel has 247.5° < glo7i < 250°. 12 contours are shown 
at logarithmic intervals with the lowest contour at a count of 5 stars per 0.1 mag x 0.025 color 
bin and the highest contour at the maximum of the plot. A streak in the contours which could be 
associated with the Perseus spiral arm is shown as a dotted line over the plots from 220-235° . b) 
Same as a) but for 0° < glat < 2°. c) Same as b) but for if ^ (mag) versus J — H color. 



Fig. 4. — This flipped B-band image of MlOO, originally published by de Jong (1996), was retrieved 
from the NED extragalactic database. Note that for an observer in MlOO looking outward away 
from the Galactic Center (as we are when we look at galactic longitudes between 220 and 250°), 
we expect to see first extinction followed by larger numbers of young blue stars as a function 
of distance from the observer. Extinction will move the luminosity function to the right in the 
color mag diagram shown in Figures 1,2 whereas a blue stellar population will move it to the left, 
countering the affect of the extinction. We expect a larger effect in the color magnitude diagram 
when we look at spiral arms that are closer to the Galactic center than the sun because we expect to 
see a younger stellar population followed by extinction as a function of distance from the observer. 



Fig. 5. — a) Similar to Figure 2a but for the region 65 < glon < 85° and 2 < glat < 3°. b) as 
in a) but for Ks{mag) vs J — H. There is a lot of extinction at glon ~ 8° which is caused by the 
Dark Nebula LDN-906 (Lynds 1962) which has an area of 15° and must be quite nearby. In this 
direction we are looking between the Carina and Perseus spiral arms. In this direction even weak 
spiral should not extend over a large range in galactic longitude, so we do not expect coherent 
features that extend from one plot to another. The large number of features seen in these plots 
is consistent with a patchy ISM and the possibility that the outer Milky Way is flocculent in its 
morphology. 



Fig. 6. — a) Similar to Figure 2a but for the region 170 < glon < 190° and —1 < glat < 1°. b) 
as in a) but for —3 < glat < —1°. Between 170 < glon < 177 we see features at Ks ~ 13.8 mag 
which would be consistent with extinction located at the expected distance of the Perseus spiral 
arm at 2 kpc from the sun. Between 175 < glon < 190 we see features at Kg ~ 14.5 which would 
correspond to a band of extinction at about 3 kpc. There is no doubt that features are coherent 
in the sense that they tend to be visible in the color magnitude plots in neighboring regions. We 
do not see evidence for one strong band of extinction associated with a dominant spiral arms but 
rather evidence for multiple weaker structures. 



